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ABSTRACT 

A millimeter molecular line survey of three carbon-rich AGB stars and two oxygen- 
rich planetary nebulae has been carried out over the frequency range 80.5-115.5 GHz. 
Sixty eight different transitions were detected in the data from 27 different molecular 
species. The hyperfine structure of G 2 H and G^^CH has been fitted to constrain the 
optical depth of their transitions. All other transitions have been constrained on the 
basis of their line profile shapes. Rotation temperatures and column densities have 
been calculated for all possible species, with adaptations to the methods applied in 
order to account for the hyperfine structure of various transitions. From the column 
densities, carbon, silicon and sulphur isotopic ratios have been determined. The results 
corroborate IRAS 15194-5115 as a J-type star, whilst excluding IRAS 15082-4808 and 
IRAS 07454-7112 as such. 

Key words: stars: AGB and post-AGB, planetary nebulae: general, stars: abun¬ 
dances, stars: carbon, line: profiles, ISM: abundances 


1 INTRODUCTION 


Asymptotic giant branch (AGB) stars and planetary neb¬ 
ulae (PNe) are rich sites of dust and molecule formation. 
As a result, their spectra are full of a variety of spectral 
lines. A vast quantity of molecular species present in these 
objects emit in the millimeter region of the spectrum. From 
these spectral lines, a host of information can be gleaned, in¬ 
cluding the abundances and temperatures of species present 
as well as their extent in the envelopes of the host stars. 
Thus millimeter molecular line surveys are an ideal method 
of identifying the chemical composition and thermodynamic 
properties of the regions in which these species exist. 

J-type stars are a sub-set of carbon-rich AGB stars with 
an unusually low ratio (< 10). It is thought that 

these could be the origin of presolar A+B grains: silicon 
carbide (SiC) grains with pre-Solar System origins, found in 


meteorites which do not have a confirmed origin (see Zinner 
|2003| for further details about the SiG presolar grains). By 
identifying and subsequently studying J-type stars in detail, 
we can ascertain the isotopic characteristics of these stars 
and compare them to the measured isotopic abundances 
found in these mysterious A+B grains. 

Isotopic abundance calculations are most easily carried 
out using optically thin rather than optically thick transi- 
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tions to measure the relative abundances of two isotopo- 
logues. Moderately optically thick lines may be used but 
would require an optical depth correction or the use of a 
radiative transfer code. However, often the method used to 
measure, for example, the optical depth of GO transitions 
requires a knowledge of the carbon isotopic ratio. Transi¬ 
tions with hyperfine structure alleviate this requirement as 
the relative intensities of the different components vary with 
increasing optical depth of the overall transition. By fitting 
these hyperfine components, the optical depth may be con¬ 
strained and thus accurate isotopic ratios determined. 


This paper details a millimeter molecular line survey 
carried out over a sample of three carbon-rich AGB stars 
and two oxygen-rich PNe. The aim of this survey is to use 
high velocity-resolution spectra to identify molecules present 
in the targets, resolve hyperfine structure where possible and 
ultimately determine reliable column densities and isotopic 
ratios for all sources observed. This line survey is intended to 
compliment those previously carried out (e.g. [Woods et al.| 


2003 Edwards fc Ziur^|2013 Ramstedt fc OlofssonpOll I 
by fully resolving the lines and subsequently analysing both 
the line profiles and hyperfine structures of the transitions 
identified. 
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Table 1. Source and coordinate list 


IRAS 

Name 

RA 

Dec 

Type 

15194-5115 

II Lup 

15 23 05.1 

-51 25 58.7 

C-AGB 

15082-4808 

RAFGL 4211 

15 11 41.9 

-48 20 01.3 

C-AGB 

07454-7112 

Al Vol 

07 45 02.8 

-71 19 43.2 

C-AGB 

14192-4355 

1C 4406 

14 22 26.3 

-44 09 04.3 

O-PN 

18021-1950 

NGC 6537 

18 05 13.1 

-19 50 34.9 

O-PN 


2 OBSERVATIONS 

2.1 Targets 

The sample comprises of five galactic sources: three carbon 
stars and two oxygen-rich planetary nebulae. Further details 
about the sources can be found in Table [T] 


2.1.1 IRAS 15194-5115 


IRAS 15194-5115, also known as II Lup, is the brightest 
carbon star in the southern hemisphere at 12 gm and the 
third brightest in both hemispheres at 12 pm, with only 
IRC-l-10216 and CIT 6 being brighter ( Nyman et al.J1993 |. 
It has a pulsation period of 575 days ( Feast et al.|2003 |. Dis¬ 
tance estimates in literature studies are between 600-1200 pc 
( Nyman &: 01ofsson|1995 Woods et al.|2003 ). The mass-loss 
rate for this star has been found from non-local thermody¬ 
namic equilibrium radiative trans fer models to be 1 x 10 ~ ^ 
M 0 yr“^ at a distance of 600 pc (Ryde et al.|l999 
et al. | 20 ro I. 


Woods 


2.1.2 IRAS 15082-4808 

IRAS 15082-4808, also known as RAFGL 4211, is also a 
mass losing, carbon-rich AGB star. Its mass loss rate has 
been derived fro m models to be approxima tely 1 x 10“® Mq 
yr“^ at 850 pc ( Groenewegen et al. [2002 1. The distance to 
this star according to the literature ranges between 640 pc 
and 1500 pc with the general consensus between 640 and 
850pc ( Woods et al.|2003 Nyman fc 01ofsson||1995 1. 


2.1.3 IRAS 07454-7112 


IRAS 07454-7112 is a mass losing, carbon-rich Mira variable 


star (Menzies et al. 20061. As a bright infrared star, this 


source has been part of several line surveys in the past as well 


as being studied in-depth by Nyman & Olofsson (19951 and 
Woods et al. (20031. The distance of this source varies in the 
literature from 710 pc to 850 pc ( Woods et al. [20(33 Menzies 
et al. 20061 , and its mass loss rate has be en modelled as 
~ 10“° Mq yr“^ (Nyman & 01ofsson|l995 I at a distance of 
1800 pc . 


2.1.4 IC 4406 

IC 4406 is an oxygen-rich, axially symmetric planetary neb¬ 
ula with a molecular envelope. It has a collimated outflow 
parallel to the major axis of the nebula, and because of this 
outflow, there is a strong dependence of the shape of the line 
profile with changing position across the nebula, as shown 


by Sahai et al. (19911. The central star temperature has 
been found to be 8 x 10^ K from three dimensional pho¬ 
toionisation modelling ( Gruenwald et al.|1997 (. This nebula 
has been shown through both observation and modelling to 
have a cylindrical structure and extends over a region ap¬ 
proximately 100” by 30” ( Faes et al.|2011 1 . 


2.1.5 NGC 6537 


NGG 6537, also known as the Red Spider Nebula, is a type 
I (rich in helium and nitrogen), bi-polar planetary nebula 
with outflows extending out to an angular size of 100 ” in the 
direction NE-SW. The central star temperature is 2 x 10® 
K. A small, opaque, dusty torus exists at 2-4” from the 
centre with an asymmetric cavity, foun d by|Matsuura et al(] 
(20051 from HST Ha and H/3 maps. Sabin et al. (20071 


show, through maps taken with the Submillimetre Gommon- 
User Bolometer Array (SCUBA), that the central region of 
the nebula is approximately 20 ” in diameter and appears to 
have a toroidal structure. 


2.2 Telescope 

The data were obtained using the 22 m Mopra telescope, sit¬ 
uated approximately 450 km from Sydney, Australia. The 3 
mm Monolithic Microwave Integrated Circuit (MMIC) re¬ 
ceiver in June 2010 was used in conjunction with the Mopra 
Spectrometer (MOPS, a digital filterbank spectrometer), in 
broadband mode for all observations. The MMIC receiver 
has a frequency range of 77-116 GHz and works in dual po¬ 
larization single-sideband mod^ 

The use of broadband mode with MOPS gives an ob¬ 
servable frequency range of 8 GHz, split into 4 overlapping 
bands, each with 2048 channels over a 2.2 GHz frequency 
range. The band edges suffer from noise, but the overlap 
given in this set up compensates well for this so little data 
is lost. Three frequency set-ups were used to observe each 
source: 84.5-92.5 GHz, 90.5-98.5 GHz and 107.5-115.5 GHz, 
with the region between 98.5 GHz and 107.5 GHz unob¬ 
served. 

The observations for this project (ID: M530) were taken 
over ten days: 18/06/2010 to 27/06/2010. The rms noise on 
a typical spectrum was 0.01 K and the velocity resolution 
was 0.91 km/s. The efficiency of Mopra at this frequency 
is ~0.5, the sensitivity is 22 Jy/K and the beam size at 90 
GHz is 38”. 


2.3 Data reduction and line identification 

The observations were made in position-switching mode. 
The on-off quotienl[^ was created from these pairs and 
their corresponding Tgys (system temperature). The result¬ 
ing spectrum is the antenna temperature corrected for atmo¬ 
spheric and systematic effects. The main beam temperature 
is obtained by dividing the corrected antenna temperature 
by the efficiency of the telescope which varies with frequency. 
The telescope efficiency is documented as 0.49 at 78 GHz and 


^ Technical summary of the Mopra radio 


narrabri.atnf.csiro.au/mopra/mopragu.pdf 
^ The quotient is defined as: T* = 


telescope: 
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0.42 at 116 GH^ Therefore, for this data, the efficiency was 
taken to be 0.49 between 84.5 GHz and 98.5 GHz and 0.42 
between 107.5 GHz and 115.5 GHz. 

The data were reduced using the ATNF Spectral Analy¬ 
sis Package (ASAP). The on/off observations were matched 
in time and the resultant spectrum formed. The two polar¬ 
isations of each set of observations were averaged, weighted 
by the system temperature. Multiple observations of each 
source in each frequency set-up were taken and thus needed 
to be coadded. However, some observations were taken dur¬ 
ing poor weather and/or suffered from severe baseline rip¬ 
pling. To remove the poor data prior to coadding, the spec¬ 
trum of each individual observation had a polynomial base¬ 
line temporarily fitted and subtracted. Those observations 
that were badly affected by the weather were flagged and 
removed from all further analysis. The remaining spectra 
(taken before the polynomial baseline subtraction) were then 
averaged in time, weighted by the system temperature. A 
polynomial baseline was htted to the averaged spectrum and 
subtracted to remove any remaining baseline rippling. The 
full spectra from all sources can be found in Appendix A of 


Smith] ( |2014'| P| 


The lines were identified using the NIST databases (Lo- 
vas||1992 Markwick-Kemper et al.1|2006l. In cases where it 


was unclear as to which species was the source of a partic¬ 
ular line, all observable transitions of each of the candidate 
species were examined for further detections. In most cases, 
this resulted in a clear identification of the line. 


3 LINE PROFILES 


The spectral lines observed from circumstellar envelopes 
(GSEs) have distinctive profiles resulting from a combina¬ 
tion of the optical depth of the emitting region and the size 
of the telescope beam compared to the size of the emitting 
region. From these shapes we can ascertain the extent of 
the emitting region relative to the beam, the motion of the 
emitting molecules and the optical depth of the emitting 
species. 

Morris ( 1975]), using the Sobolev (19601 approximation, 


derived line profile shapes for circumstellar envelopes under 
the following assumptions: 


(i) the CSE is spherically symmetric and isothermal; 

(ii) mass loss rate is constant with time; 

(iii) the expansion velocity is constant with radius; 

(iv) the local line width due to thermal motions and tur¬ 
bulence at any point in the envelope is much less than that 
created from Doppler broadening due to the expansion of 
the envelope; 

(v) the central star is negligible in size in comparison to 
the size of the envelope. 


The assumption of constant expansion velocity with ra¬ 
dius does not hold in the innermost regions of the GSE where 
the majority of the wind acceleration occurs. However, this 
region is extremely small when compared with the overall 


® Technical summary of the Mopra radio telescope: 
narrabri.atnf.csiro.au/mopra/mopragu.pdf 

https://www.escholar.manchester.ac.uk/item/?pid=uk-ac- 
man-scw:240595 


Table 2. Details of all detected transitions. All molecular tran¬ 
sitions relate to the ground vibrational state unless otherwise in¬ 
dicated. The quantum numbers of SiC 2 and C-C 3 H 2 are given in 
the assymmetric top standard format. 


Molecule 

Freq 

(GHz) 

Transition 

3°SiO 

84.74617 

J=2-l 

HCsN 

85.20160 

J=32-31 


85.22933 

N=l-0, J=3/2-1/2, Fi=2-1, F=5/2-3/2 


85.23276 

N=l-0, J=3/2-1/2, Fi=2-1, F=3/2-1/2 


85.24771 

N=l-0, J=3/2-1/2, Fi = 1-0, F=1/2-1/2 


85.25695 

N=l-0, J=3/2-1/2, Fi = 1-0, F=3/2-1/2 


85.30397 

N=l-0, 4=1/2-1/21 Fi = l-1, F=1/2-3/2 


85.30769 

N=l-0, J=l/ 2 -l/ 2 , Fi = l-1, F=3/2-3/2 


85.31438 

N=l-0, J=l/ 2 -l/ 2 , Fi=0-1, F= 1 / 2 - 1/2 

C-C 3 H 2 

85.33891 

2 i ,2 - lo,l 

C 4 H 

85.63400 

N=9-8, 4=19/2-17/2 

C 4 H 

85.67257 

N=9-8, 4=17/2-15/2 

29SiO 

85.75920 

4=2-1 

RisCN 

86.34018 

4=1-0 

SiO 

86.84700 

4=2-1 

HNi^c 

87.09086 

4=1-0 

C 2 H 

87.28416 

N=l-0, 4 = 3 / 2 - 1 / 2 , F=l-1 

C 2 H'> 

87.31693 

N=l-0, 4 = 3 / 2 - 1 / 2 , F=2-l 

C 2 H'> 

87.32862 

N=l-0, 4 = 3 / 2 - 1 / 2 , F=l-0 

C 2 H'> 

87.40200 

N=l-0, 4 = 3 / 2 - 1 / 2 , F=l-1 

C 2 H'> 

87.40717 

N=l-0, 4 = 3 / 2 - 1 / 2 , F=0-1 

C 2 H 

87.44651 

N=l-0, 4 = 3 / 2 - 1 / 2 , F=l-0 

Hi^CCCN 

88.16683 

4=10-9 

HCN 

88.63185 

4=1-0 

CCCN 

89.04559 

N=9-8, 4 = 19 / 2 - 17/2 

CCCN 

89.06436 

N=9-8, 4 = 17 / 2 - 15/2 

HCN 

89.08791 

4=1-0 F=2-l v=(0,2,0) 

HCO+ 

89.18853 

4=1-0 

HCi^CCN'’ 

90.59306 

4=10-9 

HCCi^CN'’ 

90.60179 

4=10-9 

HNC 

90.66356 

4=1-0 

SiS 

90.77156 

4=5-4 

HCCCN 

90.97899 

4=10-9 


92.49427 

4=2-1 

SiC2 

93.06364 

4o,4 - 3o,3 

SiC2 

94.24539 

42,3 - 32,2 

C 4 H 

95.15032 

N=10-9, 4 = 21 / 2 - 19/2 

C 4 H 

95.18894 

N=10-9, 4 = 19 / 2 - 17/2 

SiC2 

95.57938 

42,2 - 32,1 

C34s 

96.41295 

4=2-1 

H13CCCN 

96.98300 

4=11-10 

CS 

97.98095 

4=2-1 







unidentified 

107.973 



108.78020 

N=l-0, 4 = 3 / 2 - 1 / 2 , Fi=2-1, F=3-2 


blended transitions 
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Table 2 — continued 


Molecule 

Freq 

(GHz) 

Transition 

13CN 

108.78237 

N=l-0, J=3/2-1/2, Fi=2-1, F=2-1 

i^CN 

108.78698 

N=l-0, J=3/2-1/2, Fi=2-1, F=1-0 

i^CN 

108.78020 

N=l-0, J=3/2- 1 / 2 , Fi=2-1, F=3-2 

i^CN 

108.78237 

N=l-0, J=3/2- 1 / 2 , Fi=2-1, F=2-1 

i^CN 

108.78698 

N=l-0, J=3/2- 1 / 2 , Fi=2-1, F=1-0 

CCCN 

108.83427 

N=ll-10, 3=23/2-21/2 

CCCN 

108.85302 

N=ll-10, 3=21/2-19/2 

SiS 

108.92430 

3=6-5 

HCCCN 

109.17364 

3=12-11 

q 

00 

0 

109.78217 

3=1-0 

13CO 

110.20135 

3=1-0 

CN 

113.14419 

N=l-0, 3 = 1 / 2 -1/2) F=1/2-3/2 

CN 

113.17050 

N=l-0, 3 = 1 / 2 - 1 / 2 , F=3/2-1/2 

CN 

113.19128 

N=l-0, 3 = 1 / 2 - 1 / 2 , F=3/2-3/2 

CN*’ 

113.48812 

N=l-0, 3 = 3 / 2 - 1 / 2 , F=3/2-1/2 

CN*’ 

113.49097 

N=l-0, 3 = 3 / 2 - 1 / 2 , F=5/2-3/2 

CN*" 

113.49964 

N=l-0, 3 = 3 / 2 - 1 / 2 , F=1/2-1/2 

CN*” 

113.50891 

N=l-0, 3 = 3 / 2 - 1 / 2 , F=3/2-3/2 

CN*” 

113.52043 

N=l-0, 3 = 3 / 2 - 1 / 2 , F=1/2-3/2 

C 4 H 

114.18251 

N=12-ll, 3 = 25 / 2 - 23/2 

C 4 H 

114.22104 

N=12-ll, 3 = 23 / 2 - 21/2 

CO 

115.27120 

3=1-0 

SiCz 

115.38239 

5o,5 - 4o.4 


blended transitions 


size of the envelope. Those high energy lines formed only 
in the innermost region of the envelope will thus not fol¬ 
low the line shapes detailed in this section. Assumptions (i) 
and (ii) have been examined in the literature and clumpy 
or anisotropic structures in CSEs have been reported. This 
should be kept in mind when comparing these theoretical 
curves to observations. Assumptions (iv) and (v) are com¬ 
monly found to be true in observations. These assumptions 
do not hold for the two planetary nebula sources which are 
not spherically symmetric and have large velocity gradients. 

For optically thin emission in an unresolved source, 
the profile was shown by Morris (19751 to be flat-topped, 
whereas optically thick, unresolved emission has a parabolic 
profile. For resolved emission, these profiles become flattened 
parabola and horned profiles for optically thick and thin 
emission respectively. 


of different candidate species were identihed as potentially 
causing this transition. The Si'^^CC line at 107.971 GHz has 
been detected in IRC-l-10216 and could be the origin of this 
line, however it would have been expected that the lines 
at 97.295 GHz and 112.593 GHz which have the same and 


greater intensities respectively in IRC-l-10216 (Cernicharo 


et al. 19911. These two lines are not detected in our data. 


which leaves the identification uncertain. It is also possible 
that this transition could be an artefact of baseline ripple re¬ 
moval, although its strength and the stable baseline regions 
surrounding it make this option less likely. The line profile 
is unclear and thus constraints about the optical depth of 
this transition cannot be deduced. 

The peak intensities, integrated intensities and full line 
width were measured for every transition in all hve sources. 
For unblended transitions, the method was as follows. The 
spectra were smoothed, generally by a factor of two although 
a factor of four was used for some of the weakest lines, then 
a straight-line baseline was fitted to the data on either side 
of the line, usually using the regions offset from the line by 
±50-100 km/s, to correct for any systematic offset of the 
spectrum from zero. The limits of the line were identified 
manually, the peak intensity of the line was measured, and 
the integrated intensity was measured between the line of the 
spectrum and background. The line width was measured as 
the full width at the base of the profile. These parameters 
are noted for all transitions in all sources in Table (AGB 
star sample) and Table |^(PN sample). Absolute flux cali¬ 
bration uncertainties have been estimated at 10%. The line 
widths are accurate to 2.5 km/s. The uncertainties listed in 
the tables include both calibration and measurement uncer¬ 
tainties. 

The blended transitions fell into two categories: separa¬ 
ble blended transitions and unseparable blended transitions. 
The unseparable transitions were those that were blended 
to such an extent that the edges and thus the extent of 
the overlap of the individual lines were indeterminable. The 
integrated intensities, peak intensities and line widths for 
these blends were taken for the overall blended line and are 
indicated as such in Tables and Examples include the 
CN transitions at ~ 113.49 GHz and the HC^^CCN and 
HCG^^CN transitions at ~ 90.6 GHz, shown in Fig. Al to 
A3. 

The separable blended transitions are those that are 
on the limit of being blended, as for C 2 H, or have clearly 
defined components, as in C^^CH. These transitions have 
been fitted to measure the integrated intensity emanating 
from each component and these integrated intensities are 
listed individually in Tables and 


4 RESULTS 

Listed in Tableware details of all the transitions detected in 
the sample and a hgure for every transition in each source 
can be found in the Appendix. A total of 64 transitions were 
detected across all sources, including some blended hyperfine 
transitions. 

One emission line detected in IRAS 15194-5115 could 
not be identified. Its observed frequency is approximately 
107.973 GHz, assuming no signihcant velocity offsets with 
respect to the systemic velocity of the envelope. A number 


4.1 Line fitting 

In order to fit the line prohles of the spectra to disentan¬ 
gle the hyperhne structure of transitions and to ascertain 
the optical depth, an optimization routine is required. The 
routine chosen for this task is the Genetic Algorithm (GA) 
PiKAlA. GAs are a class of optimization routine based upon 
the principle of “survival of the fittest” employed in nature 
to produce species that adapt to their environment. For the 
initial population, the parameters are chosen using a ran¬ 
dom number generator over the entire X-dimensional pa¬ 
rameter space. In subsequent generations, the parameters 
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are constructed from the parameters of the parent solutions 
with a possibility of parameter mutation. The likelihood of a 
given individual solution being used as a parent for further 
solutions is proportional to how well it fits the data. The 
reproduction plan chosen for this work was full generational 
replacement with elitism included. Further details on this 


routine may be found in Charbonneau (19951. 



4 . 1.1 C 2 H 

C 2 H has six detected hyperfine components in IRAS 15194- 
5115 and four in IRAS 15082-4808. The profiles of C 2 H 
have an unusual appearance due to the lines being partially 
blended. The individual hyperfine components have square 
profiles, with the brightest in IRAS 15194-5115 showing a 
slight horned profile, indicating that C 2 II is just resolved in 
the beam of the telescope and is optically thin. The two- 
level step profiles are caused by two hyperfine transitions on 
the verge of blending and the three-step profiles are caused 
by two square profile transitions with a significant overlap, 
resulting in excess emission in the central region of the pro¬ 
file. 

The assumption about optical depth from the line pro¬ 
file shape can be confirmed using the method of [Fuller fc| 
Myers (19931, but in contrast to this method’s usage in that 


paper, the line profile is assumed to be a top-hat function 
rather than a Gaussian distribution. 

The optical depth of the transition, r, is now found 
using: 


t{v) = Ty^ai(f)i{v), 


( 1 ) 

where at is the relative intensity of hyperfine component i 
and (pi is the line profile of hyperfine component i, given by: 

1 if (vo + Vi — Av) < V < {vo + Vi + Av) 

0 otherwise, 

where Vi is the offset velocity of component i with respect 
to the reference velocity, vo, and Av is the line half-width. 
The line may then be fitted using: 


T(w) = A(^l-e-"'''^) +C, 


( 2 ) 


where A is related to the peak intensity of the line and C is 
a constant to fit for any residual baseline offset from zero. 
It should, however, be noted that in the optically thin limit, 
this equation reduces to T{v) = At(v) -f C, thus the ampli¬ 
tude and optical depth cannot be separately calculated in 
this limit. 

The hyperfine structure was fitted to Equ. and [fus¬ 
ing PiKAlA. The results of the fitting are shown in Fig.^ In 
both cases, the fitting suggests that the transition is opti¬ 
cally thin. The resulting integrated intensities of each of the 
hyperfine components are listed in Table 

4.1.2 

C^^CH has seven hyperfine components which, in the spec- 
trum of IRAS 15194-5115, are blended into three line- 
regions. The same approach as outlined above was applied 



Figure 1. C 2 H fitted for hyperfine structure. The upper spectrum 
is that of IRAS 15194-5115 and the lower is that of IRAS 15082- 
4808. The red line indicates the best fit to the data and the blue 
vertical lines indicate the central velocities of each of the fitted 
components. 
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Figure 2. C^^CH fitted for hyperfine structure in IRAS 15194- 
5115. 


to the C^^CH spectrum and again, the result was found to 
be optically thin. The results of the fitting are shown in Fig. 

It should be noted, however, that these lines have a low 
signal-to-noise ratio so the fits are less certain than those 
of C 2 II. The resulting integrated intensities of each of the 
hyperfine components are listed in Table 


4.2 Detections in individual sources 

4 . 2.1 IRAS 15194-5115 

IRAS 15194-5115 had the largest number of detected transi¬ 
tions of all sources in the sample. CO and both show a 
sharp peak on the blueshifted edge of the line profile, shown 
in Fig. This is caused by contamination from an interstel¬ 
lar cloud that lies in the same line of sight as IRAS 15194- 
5115, further details of which can be found in [Nyman et al.j 
(19871 and Nyman et al. (19931. When measuring the peak 


and integrated intensities of these lines, the contaminant 
emission was removed and the line interpolated between the 
remaining points. 

Transitions originating from silicon based molecules 
(e.g. SiO and SiC 2 ) have been detected, including the sil¬ 
icon isotopologues of SiO, allowing silicon isotopic ratios to 
be calculated. 

Multiple '^®C-containing species were detected in this 
source including, but not limited to, the isotopologues 
of HNC, CCH and HCCCN. Emission lines from molecules 
with a minor isotope of Si or S were also detected. 
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Table 3. Detected transitions in the AGB star sample. Frequencies are given in GHz, the line peak, T^^aki is given in K, the integrated 
intensity, f Tdv, is in K km/s and the full width of the line at the base of the line, Au, is in km/s. All temperatures are given in corrected 
antenna temperatures. The efficiency of Mopra is 0.49 before the band gap and 0.42 after. Absolute flux calibration uncertainties have 
been estimated at 10%. The line widths are accurate to 2.5 km/s. The uncertainties listed in the table in parentheses after the value 
indicate the uncertainty on the final digit and include both calibration and measurement uncertainties. 


Molecule 

Frequency 

(GHz) 

IRAS 15194-5115 

Tpeak fT'd.'IJ 

IRAS 15082-4808 
Tpeak ^Tdv 

Av 

IRAS 07454-7112 
Tpeak fTdv 

Av 

3°SiO 

84.74617 

0.014(6) 

0.3(1) 

39.1 







HCsN 

85.2016 

0.018(7) 

0.5(2) 

48.4 








85.22933 

0 . 011 ( 6 ) 

0.5(3) 

41.5 








85.23276 

0.007(6) 

0.3(2) 

41.5 








85.24771 

0.003(5) 

0 . 1 ( 2 ) 

41.5 








85.25695 

0.007(5) 

0.3(2) 

41.5 








85.30397 

0.003(5) 

0 . 1 ( 2 ) 

41.5 








85.30769 

0.006(5) 

0 . 2 ( 2 ) 

41.5 








85.31438 

0.003(5) 

0 . 1 ( 2 ) 

41.5 







C-C 3 H 2 

85.33891 

0.042(9) 

1 . 0 ( 2 ) 

43.6 







C 4 H 

85.634 

0.043(9) 

1.4(3) 

43.4 







C 4 H 

85.67257 

0.045(10) 

1.5(3) 

41.5 







29SiO 

85.75920 

0.018(7) 

0.5(2) 

40.5 








86.34018 

0.52(6) 

17(2) 

55.2 

0.07(1) 

1.8(3) 

41.2 

0.062(9) 

1.4(2) 

34.6 

SiO 

86.847 

0.17(2) 

5.3(7) 

46.5 

0.08(1) 

1.8(3) 

37.3 

0.039(7) 

0.7(1) 

27.9 

HNi^C 

87.09086 

0.024(7) 

0 . 6 ( 2 ) 

48.3 







C 2 H 

87.28416 

0.013(6) 

0.5(2) 

40.1 







C 2 H'> 

87.31693 

0.13(2) 

5.1(7) 

40.1 

0.04(1) 

1.5(4) 

36.9 




C 2 H'> 

87.32862 

0.06(1) 

2.4(4) 

40.1 

0 . 021 ( 8 ) 

0.8(3) 

36.9 




C 2 H'> 

87.402 

0.08(1) 

3.0(5) 

40.1 

0 . 021 ( 8 ) 

0.8(3) 

36.9 




C 2 H'> 

87.40717 

0.028(8) 

1.1(3) 

40.1 

0.009(7) 

0.3(2) 

36.9 




C 2 H 

87.44651 

0.014(6) 

0.5(2) 

40.1 







Hi^cCCN 

88.16683 

0.028(8) 

0.9(2) 

44.0 







HCN 

88.63185 

0.51(6) 

15(2) 

56.5 

0.44(5) 

10 ( 1 ) 

43.8 

0.14(2) 

2.8(3) 

35.6 

CCCN 

89.04559 

0.022(7) 

0.7(3) 

44.0 

0.024(8) 

0.4(1) 

37.8 




CCCN 

89.06436 

0.023(7) 

0.7(3) 

42.6 

0.024(8) 

0 . 6 ( 2 ) 

38.7 




HCN 

HCi^CCN*' 

HCCi^CN* 

89.08791 
90.59306 1 
90.60179 J 

0.043(9) 

1.8(4) 

74.0 

0.67(7) 

1 . 1 ( 1 ) 

6.3 




HNC 

90.66356 

0 . 10 ( 1 ) 

3.2(5) 

44.6 

0.05(1) 

1.4(3) 

40.1 

0.010(4) 

0.3(1) 

69.5 

SiS 

90.77156 

0.06(1) 

2.1(4) 

44.5 

0.05(1) 

1.6(4) 

45.4 

0.023(5) 

0.5(1) 

36.5 

HCCCN 

90.97899 

0 . 10 ( 1 ) 

3.2(5) 

44.4 

0 . 10 ( 1 ) 

2.7(4) 

39.1 

0.049(8) 

0.9(1) 

28.4 

13CS 

92.49427 

0.05(1) 

1.8(3) 

44.6 

0.014(7) 

0.16(8) 

34.1 




SiC 2 

93.06364 

0.07(1) 

1.9(3) 

42.5 

0.045(7) 

0.9(2) 

38.2 

0.024(5) 

0.41(9) 

31.3 

SiC 2 

94.24539 

0.040(9) 

1 . 1 ( 2 ) 

42.0 

0.026(9) 

0 . 6 ( 2 ) 

44.1 

0.014(4) 

0 . 20 ( 6 ) 

30.0 

C 4 H 

95.15032 

0.045(9) 

1.5(3) 

45.0 

0.025(9) 

0.5(2) 

41.6 




C 4 H 

95.18894 

0.047(9) 

1.6(3) 

43.3 

0.026(9) 

0.5(2) 

41.6 




SiC 2 

95.57938 

0.034(8) 

1 . 0 ( 2 ) 

43.1 

0.025(9) 

0 . 6 ( 2 ) 

38.0 

0.012(4) 

0.18(6) 

26.2 

C34S 

96.41295 

0.025(8) 

0.7(2) 

43.6 

0.012(7) 

0.3(2) 

36.9 




Hi^cCCN 

96.983 

0.017(7) 

0.4(1) 

43.3 







CS 

97.98095 

0.31(4) 

10 ( 1 ) 

42.1 0.21(2) 

- Band Gap- 

6.4(8) 

41.2 

0 . 10 ( 1 ) 

2.0(3) 

28.9 

unidentified 

iScN'^ 

107.973 

108.78020 

0.023(7) 

0.3(1) 

38.2 







i3cN*> 

i3cN*> 

108.78237 

108.78698 

> 0.027(8) 

0 . 6 ( 2 ) 

58.7 







CCCN 

108.83427 

0.018(7) 

0.5(2) 

48.3 

0.032(9) 

0.5(2) 

37.9 




CCCN 

108.85302 

0.026(8) 

0 . 8 ( 2 ) 

47.5 

0.027(9) 

0.7(2) 

36.4 





blended transitions 
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Table 3 — continued 


Molecule frequency 
(GHz) 

IRAS 15194-5115 

Tpeak 7 Tdv Av 

IRAS 15082-4808 

Tpeak 7 Tdv Av 

IRAS 07454-7112 

Tpeak 7 Tdv Av 

SiS 

108.9243 

0.045(9) 

1.3(3) 

44.5 

0.05(1) 

1 . 2 ( 2 ) 

35.6 

0.031(6) 

0.5(1) 

29.7 

HCCCN 

109.17364 

0.06(1) 

1.9(4) 

44.4 

0.09(1) 

2.5(4) 

40.7 

0.057(9) 

1 . 1 ( 2 ) 

31.8 

i^CO 

110.20135 

0.27(3) 

8.5(10) 

46.2 

0.06(1) 

1.4(3) 

40.3 

0.08(1) 

1 . 6 ( 2 ) 

29.3 

CN 

113.14419 

0 . 012 ( 6 ) 

0.8(4) 

31.4 

0.06(1) 

1 . 1 ( 2 ) 

38.5 

0.02(5) 

0.28(7) 

32.9 

CN 

113.1705 

0.05(1) 

1 . 1 ( 2 ) 

43.6 

0.06(1) 

0.9(2) 

39.3 

0.044(7) 

0 . 6 ( 1 ) 

28.6 

CN 

113.19128 

0.06(1) 

0.7(1) 

30.0 

0.07(1) 

0.9(2) 

32.9 

0.048(8) 

0 . 6 ( 1 ) 

30.7 

CN*> 

113.48812 










CN*" 

113.49097 










CN*" 

113.49964 

• 0 . 11 ( 2 ) 

4.5(7) 

112.5 

0.17(2) 

8 ( 1 ) 

127.4 

0 . 11 ( 1 ) 

3.3(4) 

77.6 

CN*" 

113.50891 










CN'> 

113.52043 










C 4 H 

114.18251 

0.06(1) 

1.5(3) 

39.6 

0.026(8) 

0 . 6 ( 2 ) 

43.2 




C 4 H 

114.22104 

0.07(1) 

1.8(3) 

42.4 

0.030(9) 

0.7(2) 

43.2 




CO 

115.2712 

0.49(5) 

19(2) 

51.2 

0.60(7) 

19(2) 

41.4 

0.52(5) 

11 ( 1 ) 

28.7 

SiC 2 

115.38239 

0.044(9) 

1 . 0 ( 2 ) 

49.0 

0.031(9) 

0 . 6 ( 2 ) 

37.8 

0.027(6) 

0.32(7) 

25.9 


blended transitions 


Table 4. Detected transitions in the PNe sample. Frequencies are given in GHz, the line peak, Tpeafei is given in K, the integrated 
intensity, f Tdv, is in K km/s and the full width of the line at the base of the line, An, is in km/s. Absolute flux calibration uncertainties 
have been estimated at 10%. The line widths are accurate to 2.5 km/s. The uncertainties listed in the table in parentheses after the 
value indicate the uncertainty on the final digit and include both calibration and measurement uncertainties. For the CO isotopologues 
in NGC 6357, the integrated intensity is that of the major peak only due to severe absorption being present in all these profiles. The 
linewidth of the CO isotopologues is across the whole profile, including the absorption region and secondary peak. 


Molecule 

Frequency 

(GHz) 

'^peak 

IC 4406 

7 Tdv 

Au 

NGC 6357 
Tpeak 7 Tdv 

Av 

Ri^CN 

86.34018 

0.025(9) 

0.4(1) 

52.4 

0 . 021 ( 8 ) 

0.3(1) 

24.3 

HCN 

88.63185 

0.04(1) 

0.9(2) 

57.4 

0.08(1) 

0.9(2) 

20.1 

HCO+ 

89.18853 

0.03(1) 

0 . 8 ( 2 ) 

55.3 

0.04(1) 

0.5(1) 

18.1 

HNC 

90.66356 

0.03(1) 

0.4(1) 

46.3 

0.06(1) 

0.43(9) 

14.3 

i^CS 

92.49427 




0 . 020 ( 8 ) 

0.18(7) 

17.5 

C34s 

96.41295 




0.027(8) 

0.25(8) 

17.6 

CS 

97.98095 




0.05(1) 

0.5(1) 

18.1 















ci®o 

109.78217 




0.07(1) 

0.35(7) 

49.3 

i^CO 

110.20135 

0.03(1) 

0.4(1) 

53.5 

0.40(5) 

2.3(3) 

45.5 

CN'’ 

113.14419 ] 









\ 0.05(1) 

0.32(7) 

24.3 




CN'’ 

113.1705 I 







CN'’ 

113.48812 1 







CN'’ 

113.49097 









\ 0.07(1) 

2.1(4) 

101.8 




CN'’ 

113.49964 







CN'’ 

113.50891 J 







CO 

115.2712 

0.16(2) 

3.2(5) 

51.9 

0.50(6) 

1.3(2) 

52.6 


blended transitions 


As explained in Sect.[^ the extent of the emission and 
the optical depth of the transitions may be constrained by 
examining line prohles. CO, for example, has a fiat topped 
parabolic profile with the counterpart displaying a 

horned profile, implying optically thick and thin resolved 
emission respectively. 

HCN displays a Gaussian line shape, implying it does 


not originate from a region of constant velocity expansion. 
This line, as shown in Smith et al. (20141, can be well mod¬ 
elled by a blend of the three hyperhne components, each 
having a Gaussian velocity distribution. 

The hyperhne structure lines in C4H are blended, with 
their relative intensities impossible to determine from the 
observed blended transition. Therefore the relative intensi- 
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Figure 3. (upper) and (lower) emission lines in IRAS 

15194-5115. Dotted lines indicate the interstellar component. Or¬ 
dinate axis intensities are in units of corrected antenna temper¬ 
ature, abscissa values are in units of km/s, corrected for LSR 
velocity of the source, taken as -15.2 km/s. 



Figure 4. HCN maser line in IRAS 15082-4808. Ordinate axis 
intensities are in units of corrected antenna temperature, abscissa 
values are in units of km/s, corrected for LSR velocity of the 
source, taken as -3 km/s from|Groenewegen et ah] ||2002[|. 



ties were calculated according to Equ. Only the major 
hyperfine components, which make up ~ 99% of the transi¬ 
tion intensity have been considered: for example, the N= 9-8, 
J=i9/2 — fins fine is considered to be split amongst the 
F=9-8 and F=10-9 components only. The blended emission 
lines are square in profile, implying the assumption that the 
transitions are optically thin is valid. However, this also im¬ 
plies that the C 4 H emission region is unresolved, although it 
is unlikely to be highly unresolved as numerous other species 
display resolved line profile shapes. 

Two transitions of HC 3 N were detected and these are 
highly blended hyperfine transitions. Einstein coefficients, 
upper level energies and statistical weights are available for 
the overall J=10-9 and J=12-ll transitions, thus alleviating 
the need to estimate the relative intensities of the unresolved 
hyperfine lines. In this source, the J=12-ll emission line has 
a square profile, suggesting optically thin, unresolved emis¬ 
sion. The J=10-9 profile is less clearly defined, with straight 
edges but a more curved central region, suggesting this line 
could be bordering on optically thick and unresolved. 


4 . 2.2 IRAS 15082-4808 

IRAS 15082-4808 had the second highest number of tran¬ 
sitions of the sample. The most unusual detection made 


Figure 5. CN emission lines in IRAS 15082-4808. Ordinate axis 
intensities are in units of corrected antenna temperature, abscissa 
values are in units of km/s, corrected for LSR velocity of the 
source, taken as -3 km/s from|Groenewegen et ah] ||2002[|. 


is that of the vibrational (v 2 = 2 ) HCN maser emission at 
89.087 GHz. Only ten of these masers, including this de- 


tection, have been found to date (Smith et al.| 

2014 

Lucas 

et al.|1986| |Guilloteau et al.||1987| |Lucas et al. 

1988 

Lucas 


found in|Smith et al. (20141 and the maser line is shown in 

Fig. El 

The CN line at 113.499 GHz shows interesting struc¬ 
ture. The blue shifted part of the spectrum shows a signif¬ 
icantly shallower drop than that on the red. This is due to 
a blend of hyperfine lines causing the appearance of a shal¬ 
lower drop. Similarly the HC 3 N 109.173 GHz line has this 
feature, but on the red shifted region of the spectrum and is 
also likely to be caused by a hyperfine blend. The hyperfine 
structure lines in C 4 H are blended and treated in the same 
fashion as in IRAS 15194-5115. The C 4 H lines have a low 
signal-to-noise ratio which makes it difficult to classify the 
line profiles. They appear to have noisy square profiles, but 
this is not certain. 

The line profile of these HC 3 N transitions in IRAS 
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15082-4808 are parabolic, implying that these lines are op¬ 
tically thick and unresolved. 

Isotopologues of a number of species were detected in 
this source, including the and isotopologues of CS, 
CO and HCN and the and isotopologues of CS. 


4.2.3 IRAS 07454-7112 

IRAS 07454-7112 has the lowest number of detected tran¬ 
sitions in the AGB star portion of the sample. The lines 
have significantly lower intensities than in the previous two 
AGB stars. Many of the detected transitions have parabolic 
shapes and the remainder have too low a signal-to-noise ra¬ 
tio to reliably ascertain the shape of the profile. Despite 
having similar Tpeak for GO, HGN and CN lines to the 
aforementioned stars, C 3 N and the polyyne radicals G 2 H 
and G 4 H were not detected down to an RMS noise level 
of 0.01 K. However, three silicon bearing molecules were 
detected, including three transitions of SiG 2 . Similarly to 
IRAS 15082-4808, two transitions of HG 3 N were detected 
as blended transitions, each with parabolic profiles implying 
optically thin emission.. 


4.2.4 ic 4406 


IC 4406 has the lowest number of detections of all of the sam¬ 
ple, with transitions from only seven species detected. This 
is to be expected as planetary nebulae are less molecule-rich 
than AGB stars due to the high levels of ionising radiation. 

CO has a double-peaked appearance in this target; the 
blue-shifted peak is more intense than the red-shifted peak. 
This nebula is not spherically symmetric, thus it is difficult 
to infer information about optical depth and whether or not 
the nebula is resolved from the line profiles. The variation in 
line prohles of CO with position on the nebula were mapped 


by Sahai et al. (19911 and this shows the positional depen¬ 


dence on the profile shape, caused by different gas motions. 
As the double-peaked appearance and relative intensities of 
the peaks are strongly dependent on observation position on 
the nebula, it appears that the double-peaked appearance is 
likely due to morphology. The ^^CO line shows only the blue- 
shifted component, explained by its significantly less intense 
lines and the exact position of the telescope on the source. 
The integrated intensity of CO in Table is the sum of the 
integrated intensities of the two components. 

The remaining transitions also show the double-peaked 
appearance, although in the isotopologues of HCN and in 
HCO^, the red-shifted component is more intense than the 
blue shifted peak. 


4.2.5 NGC 6537 


5 ANALYSIS 


The analysis carried out on the data consists of: the mea¬ 
surement of the rotational temperatures of different species 
using the population diagram method, the calculation of 
source-averaged column densities of each species and iso¬ 
topic ratios from a number of sets of isotopologues. For the 
following analysis, the antenna temperature has been con¬ 
verted to main beam temperature, assuming an efficiency 
of 0.49 before the band-gap and 0.42 afterwards. All the 
constants used in the analysis originate from the JPL and 
CDMS catalogues ( Pickett et al.|[r998 Muller et al.||2005 |. 
Unless otherwise stated, lines are assumed to be optically 
thin. 


5.1 Excitation temperatures 


For those species where multiple transitions were detected 
(excluding hyperfine transitions), it was possible to use the 
population diagram method, as presented in [Goldsmith fc| 
Langer (|1999 1, to measure the rotational temperatures of 


different species. It has been assumed that the excitation 
temperature can be taken as the rotational temperature in 
the following analysis. The details of the method used are 
outlined below, including adaptations to account for hyper¬ 
fine structure. 


5.1.1 Standard population diagram method 

The population of the upper level of a transition, Nu, can 
be expressed as: 


= (3) 

where N is the column density of the molecule, Eu is the 
energy of the upper level, Qu is the statistical weight of the 
upper level and Tex is the excitation temperature. The par¬ 
tition function, Q is given by: 


« = E^* = E5»®-p(-i^)- (4) 

where the sum is over all energy levels, i. 

Thus, rearranging and taking logs of Equ. [pleads to: 

which, when plotting In against gives a straight 

line with gradient equal to —jt— and y intercept equal to 
(In A — InQ). The population of the upper level, assuming 
the source fills the beam and the transition is optically thin, 
is given by: 


NGC 6537 has transitions detected from 10 species, although 
all three CO isotopologues (^^CO, ^®CO and C^®0) show 
probable interstellar contamination - these line regions can 
be seen in Fig. A5. 

Two isotopologues of HCN were seen (HCN and 
H^^CN) along with HNC. Three isotopologues of CS were 
detected: CS, ^^CS and C®'*S. HCO^ was also detected. 


_ 87rkBi^^IU 
“ ' hc3A,i 


( 6 ) 


where W is the integrated intensity of the line and Aui is 
the Einstein A coefficient of transition u-l. 

For some species with hyperfine structure, the input 
parameters (e.g. A„i) for the population diagram method 
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are available for the overall N or J transition (e.g. HCCCN). 
In these cases, Equ.j^may be directly applied, taking W as 
the total integrated intensity of all hyperfine transitions. 

In many cases, the input parameters are only available 
for the individual hyperfine components. In this case, the 
method outlined above must be adapted and these adapta¬ 
tions are detailed in the following subsections. 


5.1.2 Hyperfine structure: separable 


If all hyperfine components of the overall N or J transition 
have been detected and they are either separated or sepa¬ 
rable (as is the case for C 2 H in IRAS 15194-5115, see Sect. 


4.1.1 1 so the integrated intensities of all individual hyperfine 


components are known, the population of the upper level can 
be taken as the sum over all components: 


Nu 


SrrkB 

hc3 


E 



(7) 


where the summation is over all hyperfine transitions, i. 

If all components are not detected, but those that are 
detected are separated or separable, the method of calculat¬ 
ing the population of the upper level must be adapted to 
compensate for those components not seen. Thus the upper 
level population becomes: 


N„, = 


87rkB 1 
hc^ y ] Uj^obs ' 


E 


i,obs 


^i,obs^i,obs 

^i,obs 


( 8 ) 


with the summations running over all observed hyperfine 
transitions (i, obs) and ai,obs are the relative intensities of 
the observed hyperfine transitions. For this work, the rel¬ 
ative intensities of the lines have been taken as the rela¬ 
tive intrinsic line strengths, taken from the JPL or CDMS 
databases and normalised ( [Pickett et al.|[T998| [Muller et al.j 
2005[). Th ese parameters can also be found in Appendix B 
of Smith ( [2014[ p| 

It is important to note that there is a distinct difference 
between the upper level, u, which is the unsplit upper energy 
level, and that of the individual hyperfine transitions, i. In 
cases where En, the energy of the unsplit upper level, is not 
known, the mean of the energies of the upper levels of the 
hyperfine transitions, ^ jj^s been taken in its place. 


5.1.3 Hyperfine structure: unseparable 

For those transitions that are blended to the extent that they 
cannot be disentangled, yet the required input parameters 
for the population diagram method are only available for 
the individual components, a further adaptation is required. 
It has been assumed that the integrated intensities of the 
blended transitions are divided amongst their components 
according to the intrinsic line strength: 


® https://www.escholar.manchester.ac.uk/item/?pid=uk-ac- 
man-scw:240595 


Table 5. Calculated excitation temperatures for a variety of 
species (K). 


Species 

IRAS 

15194-5115 

IRAS 

15082-4808 

IRAS 

07474-7112 

C 4 H 

79 ±56 

- 82 ± 147 


HC 3 N 

14 ±5 

35± 28 

240 ± 1260 

Hi^CCCN 

5 ± 3 



SiS 

7± 3 

10 ±6 

55 ± 288 

SiC 2 

19± 11 

35± 51 

14± 16 

C 3 N 

31 ± 21 

543 ± 5760 



m = IFblend^^. (9) 

E“^' 

i=o 

where Wi is the integrated intensity of hyperfine component 
i, ITbiend is the integrated intensity of the blend of n hy¬ 
perfine transitions including hyperfine transition i and Ui is 
the relative intrinsic line strength of hyperfine component 

71 

i. Uj is the sum of the intrinsic line strengths of all hy- 
j=o 

perfine components that make up the blend. The integrated 
intensities of the individual transitions may now be used in 
either Fqu.[7]or|^for those where all hyperfine components 
or some hyperfine components are detected, respectively. 


5 . 1.4 Line profiles and excitation temperature results 

The calculated excitation temperatures from the population 
diagram method are given in Table Population diagrams 
are shown for IRAS 15194-5115, IRAS 15082-4808 and IRAS 
07454-7112 in Fig.|^ 

Additionally, as only two transitions of this molecule are 
observed, the reliability of this method is also significantly 
decreased. The calculated excitation temperatures are given 
in Table [5] 

Two transitions of SiS (J=5-4 and J=6-5) were detected 
in all three of the AGB stars in the sample. These lines have 
no fine or hyperfine splitting associated with them. The two 
lines have a relatively square profile shape in IRAS 15194- 
5115, although they show a degree of asymmetry, implying 
that these transitions are optically thin. In IRAS 15082- 
4808, the J=6-5 line has a square profile and the J=5-4 
transition has a less clear line profile, with resemblance to 
both square and parabolic, implying these are likely to be 
optically thin. Both lines of SiS have parabolic profiles in 
IRAS 07454-7112, implying that both transitions are opti¬ 
cally thick. 

Up to four SiC 2 transition lines were detected in the 
spectra of each of the AGB stars in the sample. The lines 
in IRAS 15194-5115 and IRAS 15082-4808 appear to all be 
square, with only the ~ 95 GHz transition in IRAS 15082- 
4808 potentially being parabolic. The lines in the spectra of 
IRAS 07454-7112 have too low a signal-to-noise ratio to be 
clearly classified. 

G 3 N was detected in two objects in the sample (IRAS 
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Figure 6. Population diagrams for IRAS 15194 — 5115 (upper 
plot), IRAS 15082 — 4808 (middle plot) and IRAS 07454 — 7112 
(lower plot). All molecules are shown and identidified by their 
symbols and colours. The black dashed line indicates the fit to all 
data. 


15194-5115 and IRAS 15082-4808) with two transitions. In 
both cases, the observed line profiles are square. As with 
C 4 H, the lines are blends of unseparable hyperfine compo¬ 
nents with major components and essentially negligible mi¬ 
nor components. Thus the same method for integrated inten¬ 
sity determination has been applied to C 3 N as was applied 
to C 4 H. 


5.2.1 Column densities: population diagram method 

Column densities can be calculated from the y-intercept, C, 
of the population diagram (Equ. [^: 


C = (InAE-lnQ), (10) 

where N is the source-averaged column density and Q is 
the partition function of the molecule. The partition func¬ 
tions have been taken from the Jet Propulsion Laboratory 
(JPL) database or the Cologne Database for Molecular Spec¬ 
troscopy (CDMS) via Splatalogue ( [Remijan et al.||2007[ ) 
and are calculated for a variety of excitation temperatures. 
The partition function evaluated at the closest tempera¬ 
ture to that found from the population diagram analysis 
was adopted for each molecule, barring C 4 H in IRAS 15082- 
4808 for which a temperature of 37.5 K was assumed as the 
measured excitation temperature was unphysical. The val¬ 
ues were not interpolated as they are generally within the 
measured uncertainties of the excitation temperatures. 

The calculated column densities for all molecules for 
which population diagrams could be constructed are shown 
in Table m 


5.2.2 Column densities: single-line method 

Population diagrams, as described above, could not be cre¬ 
ated for all molecules, thus the column densities were calcu¬ 
lated for the remaining molecules by combining Equ. and 

El 


87rkBi2"VE Q f Eu 

hc^Aui gu V^BTej 


( 11 ) 


where N is the column density of the species, o is the fre¬ 
quency of transition, W is the integrated intensity of the 
emission line, Aui is the Einstein coefficient of the transition 
from the upper level, u, to the lower level, 1. Blended lines 
and hyperfine transitions were treated in the same fashion 


as Sect. 5.1 similar to that of Greaves & White (19921. If 


the parameters are available for the overall transition, Equ. 
[m becomes simply: 


N = 


87rkBi2^ Q 

hc^Auj gu 



( 12 ) 


If parameters are only available for the individual hyperfine 
transitions, Equ. |11| becomes: 


_ 87rkB Q f Eu \ V G 


w. 


(13) 


if all hyperfine components are observed, or if only some 
hyperfine components are observed: 


5.2 Column densities 

Two different methods of calculating the column densities 
of molecules have been applied. These both use the same 
assumptions: the source fills the beam, the emission is opti¬ 
cally thin and LTE applies. 


^ _ 87rkB Q f ^ 1 ^i,obsll^i,oba 

hc^ Qu VksTex / Ui obs , Ai obs 

^ — i,ohs 

(14) 

The summations run over either all hyperfine transitions, i, 
or all observed hyperfine transitions (i,obs). 

The excitation temperature was assumed as: 18.75 K 
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for IRAS 15194-5115 and 37.5 K for all other sources. These 
were selected as closest to the mean of the calculated exci¬ 
tation temperatures for each of the AGB stars (excluding 
outliers) and were simply assumed for the PN as no excita¬ 
tion temperatures could be calculated. 

The aforementioned excitation temperatures are given 
to a higher degree of precision than is realistic for our re¬ 
sults from the population diagram analysis. The tempera¬ 
tures have been specihed at that level of precision solely 
due to the fact that the partition functions, available from 
CDMS, have been calculated at these very specific tempera¬ 
tures. For the effect of the excitation temperature uncertain¬ 
ties on the isotopic ratios and column density calculations, 
see Sect. 16.71 

The results of these calculations are shown in Table ID 


0.10 

0.08 

„ 0.06 
£• 0.04 

S 0.02 
^ 0.00 
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5.3 Isotopic ratios 

Isotopic ratios are most directly measured using observa¬ 
tions of optically thin lines. The optical depth can be some¬ 
what constrained by examining the line profiles and making 
the assumptions described in Sect. 

The isotopic ratios have been calculated by comparing 
the calculated source-averaged column densities. The results 
are shown in Tabled 

5.3.1 Carbon isotopic ratios 

A number of carbon-bearing species and their correspond¬ 
ing isotopologues were identified in these observations: CO, 
HCN, HNC, CS, HCsN, CN and all of their corresponding 

isotopologues. Additionally, CCH and one of its car¬ 
bon isotopologues, C^^CH was identified. The transitions of 
the second isotopologue, ^^CCH, are outside the frequency 
range of these observations. 

For those species with a single isotopologue, such as 
CO and CS, the isotopic ratio is found by comparing the 
column densities of the two isotopologues. For those with 
multiple isotopologues, the isotopic ratio can be determined 
from each pair and should be equal, assuming there is no 
chemical fractionation towards one isotopologue. For exam¬ 
ple, HCsN has three isotopologues (H^^CCCN, HC^^CCN 
and HCC'^^CN) and the isotopic ratio may be found by di¬ 
viding the column density of the ^^C-containing isotopo¬ 
logue by each of the column densities of the ^®C-containing 
isotopologues resulting in three ^^C/^^C isotopic ratios. 

CO and HCN appear to be optically thick in all the 
ACB stars in the sample, thus the ratios found from these 
molecules are lower limits only and, as expected, give the 
lowest calculated ratios. HNC, HC 3 N and C 2 H appear to be 
optically thin from the appearance of their line profiles in 
IRAS 15194-5115 and their ratios agree well with one an¬ 
other. ^^CS in some cases appeared to be on the borderline 
between optically thick and optically thin. In IRAS 15194- 
5115, the measured isotopic ratio is lower from CS than from 
HNC and C 2 H, which would be expected if CS were some¬ 
what optically thick. The hyperfine transitions of CN were, 
in general, blended making it difficult to identify whether 
the transitions were optically thick or thin. Fig. [T] shows the 
hyperhne structure of CN with examples of both optically 
thick and optically thin hts to the data. The optically thick 


0.10 
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0.04 
J 0.02 
^ 0.00 


-200 -100 0 100 
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Figure 7. Figures showing example fits to the CN emission in 
IRAS 15194-5115. The upper two plots are optically thick fits to 
the data and the lower two are optically thin fits to the data. 

blended profile better represents the data, suggesting that 
the transition is optically thick. The calculated carbon iso¬ 
topic ratio is lower than the ratios found from HNC and C 2 H 
which also suggests that the transition is optically thick. 

The ratios calculated here are lower than those found 
by directly evaluating the integrated intensities of the lines, 
without compensating for the transitions’ different Einstein 
coefficients. For IRAS 15194-5115, for example, the inte¬ 
grated intensity ratios of isotopologues HNC and CS are 
5.5 and 5.8 respectively. 


5.3.2 Sulphur isotopic ratios 

Three isotopologues of CS were detected in multiple sources: 
i2c32s, i 2 Q 34 g i3Q32g ^ further two isotopologues 
were not detected in any of the sample stars: '^^C^^S and 
12c36s 

Deriving ®^S/®'’S can be done in two ways with the de¬ 
tected transitions: by comparing the integrated intensities 
of the ^^C^^S and ^^C^'^S lines or by dividing the integrated 
intensity of ^^C^^S with that of ^^C^^S and multiplying by 
the carbon isotopic ratio of the source. 

I 2 c 34 g^ where detected, has an optically thin profile. 
12 c32s, however, showed a borderline optically thick/thin 
profile. In IRAS 15082-4808, both methods give consistent 
ratios of approximately 22 for ^^S/^'^S. In IRAS 15194-5115, 
the ratios differ by a factor of ~ 1.5, suggesting that CS is 
indeed optically thick in this star. 


aAAi 
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Table 6. Column densities of all detected species in all sources. Quoted values are in cm“^. Uncertainties quoted here are the formal 
errors propagated from calibration and measurement uncertainties. A discussion on the effects of excitation temperature and other 
assumptions can be found in the text. Those values calculated from the rotation diagram method are shown in bold, the remaining values 
used the single-line method. The lower limits in this table originate from transitions which are likely optically thick.. 


Species 

IRAS 15194-5115 

IRAS 15082-4808 

IRAS 07454-7112 

IC 4406 

NGC 6357 

CO 

> (5.1 ± 0.5) X IQi® 

> (8.6 ±0.9) X IQi® 

> (4.8 ±0.5) X IQi® 

(1.5 ±0.3) X lOi® 



(5.0 ± 0.6) X IQi® 

(1.2 ±0.3) X IQi® 

(1.4 ± 0.2) X IQi® 

(4.4 ± 1.0) X IO15 


HCN 

> (2.2 ±0.3) X IQi'* 

> (2.7 ± 0.3) X IQi'* 

> (7.1 ± 0.8) X IO13 

(2.3 ±0.5) X IO13 

(2.3 ±0.5) X IO13 


> (2.7 ± 0.3) X IQi'* 

> (5.0 ± 0.8) X IO13 

(3.7 ± 0.5) X IO13 

(1.1 ± 0.3) X IO13 

(9± 3) X 1013 

HNC 

(1.5 ±0.2) X I013 

(1.1 ± 0.2) X 1013 

(2.6 ±0.8) X 1013 

(3.5 ±0.8) X 1013 

(3.4 ± 0.7) X 1013 

HNi^c 

(4± 1) X 10^2 





C2H 

(4.1 ± 0.4) X IQi® 

(1.5 ±0.3) X IQi® 




C^^CH 

(1.1 ± 0.4) X IQi® 





C4H 

(2.8 ± 0.7) X 10^® 

(2 ± 2) X 1014 




C3N 

(1.4 ±0.8) X 10^4 

(8 ±4) X 1014 




HC3N 

> (2 ± 1) X 10^^ 

> (3 ± 2) X 10i3 

> (3 ± 2) X 10i3 



Hi^CCCN 

{8± 20) X 10^^ 





HCi^CCN 

(1.3 ±0.3) X IO13 





HCCi^CN 

(1.3 ±0.3) X IO13 





HC5N 

(4± 1) X 10^3 





HCO+ 




(4± 1) X 1013 

(2.3 ±0.5) X 1013 

CN 

(1.2 ±0.1) X IQi® 

(2.9 ±0.2) X IQi® 

(1.4 ± 0.1) X IO13 

(6.3 ±0.6) X 10i4 


i^CN 

(4.3 ±0.9) X 10i'‘ 





CS 

> (1.1 ± 0.1) X IQi'^ 

> (1.1 ± 0.1) X 10i4 

> (3.6 ±0.5 X 1013 


(9± 2) X 10i3 

i^CS 

(4.2 ±0.7) X IO13 

(6 ±3) X 1013 



< (7± 3) X 1013 

C34S 

(8± 2) X IO12 

(5 ±3) X 1013 



(5± 1) X 1013 

C3H2 

(2.3 ±0.5) X IO13 





SiO 

> (2.7 ± 0.4) X IO13 

> (1.6 ±0.3) X IO13 

> (6.2 ±0.9) X 10i3 



29SiO 

(2.4 ± 0.9) X I0I2 





30SiO 

(1.6 ±0.7) X I0I2 





SiS 

(7± 5) X lO^^ 

(3 ± 3) X 1013 

(1 ± 2) X 1013 



SiC2 

(4 ±2) X 10^^ 

(4 ±3) X 1013 

(4 ±5) X 1013 




All three CS isotopologues were detected in NGC 6537. 
32g/34g measured to be 2 and 4 from the two aforemen¬ 
tioned methods, which is significantly lower than those of 
the AGB stars. However, the minor isotopologues had low 
signal-to-noise ratios and therefore these values should be 
taken with caution. 


5.3.3 Silicon isotopic ratios 

^®SiO, ^®SiO and ®°SiO have been detected in IRAS 15194- 
5115. ^®SiO has a clear parabolic profile, meaning only lower 
limits can be derived for ^®Si ratios. ^®SiO and ^°SiO showed 
square profiles and are thus optically thin. The measured 
29gi/30si ( 1 . 5 ) is consistent with the Solar value. 


6 DISCUSSION 
6.1 Emission extent 

From the observed line profiles, limits can be placed on the 
extent of the emitting regions. For example, the line profile 


of CO in IRAS 15194-5115 has a flattened parabolic pro¬ 
file (neglecting the contaminant emission) and has a 

double-peaked profile, both of which imply resolved emis¬ 
sion. Assuming a distance of 1 kpc to IRAS 15194-5115, 
the beam (33") resolves an area with radius 3.0 x lO^^cm, 
putting a lower limit on the CO emission region. The limit 
is dependent on the distance to the source used. Assuming a 
distance of 600 pc, the radius of the area subtended by the 
beam reduces to 1.6 x lO^^cm and if a distance of 1.2 kpc 
is assumed, this radius increases to 3.2 x 10^^ cm. The CO 
photodissociation radius reported in [Woods et al.| ( |2003| ) is 
3.2 X 10^^ cm assuming a distance of 600 pc. Our lower limit 
value is a factor of two smaller than this value but as it is a 
lower limit, it does not disagree with the results of jWood^ 


et al. (20031. 


A number of other molecules, including HC 3 N, SiO and 
their isotopologues, show unresolved profiles, implying that 
their emission region is smaller than that of CO. The up¬ 
per level energies of the resolved transitions (molecules: CO, 
^^CO, C 2 H and likely CN and CS) are less than 5 K, whereas 
the upper level energies of the unresolved emission lines are, 
as expected, higher. 
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Table 7. Isotopic ratios derived from calculated column densities. Uncertainties quoted here are the formal errors propagated from 
calibration and measurement uncertainties. A discussion on the effects of excitation temperature and other assumptions can be found in 
the text. 


Ratio 

Species 

IRAS 15194-5115 

IRAS 15082-4808 

IRAS 07454-7112 

IC 4406 

NGC 6537 

12C/13C 

CO, 13CO 

> 1.1 ± 0.2 

> 7± 2 

> 3.4 ±0.6 

3.4 ± 1.0 



HCN, Hi^CN 

> 0.8 ±0.1 

> 5± 1 

> 1.9 ±0.3 

2.1 ± 0.7 

2.6 ± 1.0 


HNC, HNi^C 

3.8 ± 1.1 






cs, 

> 2.6 ±0.5 

18 ±9 



> 1.3 ±0.6 


C2H, Ci^CH 

3.7 ± 1.4 






HC3N, H13CCCN 

> 0.3 ±0.6 






HC3N, HC13CCN 

> 1.5 ± 0.8 






HC3N, HCC13CN 

> 1.5 ± 0.8 






CN, 

> 2.8 ±0.6 





32g/34s 

c^^s, I3cs 

21 ± 6“ 

22 ± 17*’ 



4± 2= 


cs, c^^s 

> 13.8 ±3.7 

> 22± 13 



1.8 ±0.5 

28Si/29Si 

SiO, 29SiO 

> 11 ± 5 





28Si/30Si 

SiO, 30siO 

> 17±8 





29Si/30Si 

29SiO, 30SiO 

1.5 ± 0.7 







“: I2c/i3c = 4.0, *>: 

12C/13C = 18.2, 

12C/13C =2.5 




6.2 Excitation temperatures 


Rotation temperatures were determined for up to six 
molecules in the AGB star portion of the sample. IRAS 
15194-5115 has six measured rotation temper atures. These 
range from 5 K for H^^CCCN to 77 K for C 4 H. Woods et al. 
| |2003[ ) determine rotation temperatures for this source as 
part of their molecular line survey using data obtained on 
the Swedish-ESO Submillimeter Telescope (SEST) and the 
Onsala Space Observatory (OSO) 20 m telescope. Of the 
eight molecules that they were able to measure rotation tem¬ 
peratures for, three also appear in our results: H^^CCCN, 
HC 3 N and SiS. Our results agree to better than 3 K (21%) 
for HC 3 N and agree to within 1 K (15%) for SiS and ~2 
K (40%) for H^^CCCN which are all within our calculated 
uncertainties. Our calculated temperature for SiC 2 agrees 
with the average result obtained for their AGB star sample 
(24.4 K). The determined temperature of C 3 N (31 K) has 
reasonable agreement with SiC 2 . C 4 H, however, is a little 
higher at 77 K. 


Woods et al. (2003 \ have also measured rotation tem¬ 


peratures in the remaining two sources of our AGB sample. 
Our obtained value for the rotation temperatures of SiS and 
SiG 2 agree to within 5 K (50% of our value) and 2 K ( 6 % 
of our value) in IRAS 15082-4808 which are both within 
the measured uncertainties. However our calculated value 
for IRAS 07454-7112 is a factor of 10 higher than that mea¬ 
sured in Woods et al. (20031, although the uncertainties on 


our SiS excitation temperature in IRAS 07454-7112 are ex¬ 
tremely large (±288 K). 

Differences in the excitation temperature values could 
be caused by the different methods used to obtain the rota¬ 
tion temperatures or could be an effect of the different beam 
sizes of the telescopes used to obtain the data: Mopra has a 
beam size of 36”±3 at 86 GHz while SEST and OSO 20 m 
had 57 ” and 44" respectively at the time of observation. 

HC 3 N in IRAS 15082-4808 and both HC 3 N and SiS in 
IRAS 07454-7112 appear to be optically thick. This results 


in the underestimation of the upper level populations and 
will also have an effect on the rotation temperature derived 
if the optical depth of the transitions varies with upper level 
energy. 

The calculated excitation temperature of C 4 H is un¬ 
physical (negative) in IRAS 15082-4808. This is caused by 
the shallow gradient of the population diagram constrained 
by level energies that are in relatively close proximity to one 
another. This temperature has therefore not been used in 
subsequent analyses. 


6.3 Source-averaged column densities 


Source-averaged column densities have been calculated in 
each source for up to 27 different molecules. Two methods 
were used: the population diagram method and a direct cal¬ 
culation for those molecules for which population diagrams 
were not able to be constructed. The AGB stars in the sam¬ 
ple have higher numbers of molecular species detected on 
average than those of the planetary nebulae in the sample. 
This is typical for PNe in comparison to AGB stars due to 
the high levels of ionising radiation which dissociate many 
molecular species in the nebulae. 

In the work of Ali (20061, a number of molecular col¬ 


umn densities are predicted from gas phase models of IRAS 
15194-5115. Our calculated abundances for CN, HNC, SiC 2 , 
G 3 H 2 and HC 5 N agree to within a factor of two of those pre¬ 


dicted by the models of Ali (2006 1 . The abundance limits of 


molecules such as GS and HCN, whose transitions are opti¬ 


cally thick, agree with the models of Ali (20061. Our mea¬ 


sured abundance of SiS is two orders of magnitude lower 
than the predicted value. The abundance of HG 3 N is a fac¬ 


tor of 100 higher than the predictions of Ali (20061, but 


more in-line with the measured column densities of HC 3 N 
in IRC±10216 by [Woods et al.| ( |2003[ ). 

In general, the column densities of the AGB star sample 
are broadly similar, with differences of up to a factor of 3 
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between sources. Notably, IRAS 07454-7112 had no detec¬ 
tions of C 2 H, C 4 H and C 3 N, despite the column densities 
of the detected species being very similar to those of IRAS 
15082-4808 and IRAS 15194-5115. 


Edwards & Ziurys (20131 have studied NGC 6537 in- 


depth in similar frequency ranges to those observed here 
and use the non-LTE RADEX code to analyse the data. 
Our column densities for HNC agree to within 20% of that 
found by Edwards & Ziurys (20131. However, our HCN and 
H^^CN are higher than their calculated column densities by 
a factor of 4, although the ratio of the two agree to within 
4%. As the CO isotopologue regions of the spectra included 
contaminant absorption, it was impossible to calculate col¬ 
umn densities for these molecules. The work of lEdwards fcl 


Ziurys (20131 does not mention suffering from contamina¬ 


tion, although the data in this region is not shown. 


Cox et al. (19921 have studied abundances of IC 4406 


relative to HCN. Assuming all species are present over the 
same region, we can compare our source-averaged column 
densities to those abundances reported in Cox et al. (19921. 
Our relative abundance of HNC to HCN is approximately 
consistent with that found in Cox et al. (19921. However 
the relative abundances of CO, CN and HCO^ differ from 
those measured by Cox et al. (1992) by factors of 2, 1.5 and 
3 respectively. This may be due to the different extents of 
the emission regions or due to the observations having been 
taken from a different region of the nebula as the strength 
and profile of the lines has been shown to vary across the 
nebula. 

Although both nebulae are O-rich, at least in the ionised 


gas, with C/O ratios of 0.79 (IC 4406, Delgado-Inglada & 
Rodrfguez|2014 1 and 0.95 (NCC 6537, Pottasch et al.|2000 1 , 


the only oxygen containing species detected was HCO-I-. 
Species such as SO, SO 2 and SiO which are seen towards 
O-rich sources were covered by our survey, but not detected. 
The absence of SiO is perhaps not surprising as it is often 
‘locked away’ in grains in the outer regions of planetary neb¬ 
ulae (e.g. Thronson fc Lada|1982 1 . CS, an unusual species to 
be present in an O-rich system, was detected towards NGC 
6537. As mentioned in Woods et al. (20051, CS can form 
in regions of shocked gas so this may partially explain the 
presence of such a molecule. Also, as the measured C/O ra¬ 
tio is close to 1 , it is possible that at least some regions of 
this nebula could be carbon-rich. 


6.4 Isotopic Ratios 

A number of isotopic ratios have been calculated for the 
sample, assuming each of the molecular isotopologue col¬ 
umn density ratios is eqnal to the overall elemental isotopic 
ratio. Only lower limits could be applied to IRAS 07454-7112 
due to the lack of detected ^^C species, suggesting that this 
source has the highest of the sample. Taking the 

H^^CCCN line (J=10-9) as an example, the rms noise in 
the region of this line is 0.004 K (corrected antenna tem¬ 
perature units) which, assuming the line is optically thin 
and therefore square, gives an integrated intensity of 0.113 
K km/s assuming a linewidth equal to the HCCCN J=10-9 
line (28.4 km/s). This results in an upper limit on the col¬ 
umn density of 1 x 10 ^^ cm“^ and a limit on the 
of 38. 

IRAS 15082-4808, however, has been found to have a 


12c/13c j-atio of 18.2 from CS, and significantly lower limits 
from HCN and CO isotopologues. This puts this carbon-rich 
AGB star above the threshold for a J-type star. 


Both of the observed planetary nebulae are oxygen-rich. 
Their ^^C/^®C have been measured in this work at 2.4 and 
2.5 from HCN isotopologues in IC 4406 and NGC 6537 re¬ 
spectively. The value for NGC 6537 agrees to within 4% of 
that reported by [Edwards fc Ziurys| ( |2013| ) . These isotopic 
ratios are lower than those of other planetary nebulae (e.g. 


Balser et aL|2002 

1 as well as those of the oxygen-rich AGB 

stars reported by 

Milam et al.||2009|. 


IRAS 15194-5115 had the largest number of isotopo¬ 
logues detected and thus had the largest number of isotopic 
ratios calculated. The ^^C/^^C ratios are consistent across 
HNC, C 2 H and CS at 3-4. The value derived from HC 3 N and 
its isotopologues is lower than the CNO cycle equilibrium 
value, however these are lower limits. This alters, however, 
if the direct method of column density calculation is used: 
the ^^C/^^C ratio increases to ~ 4. This is the only molecule 
where a significant diference is found between using the two 
methods and it is unclear as to why this occurs. The limits 
ascertained for CO, HCN are lower than these values, thus 
are also consistent. The ^^C/^^C as derived from CN iso¬ 
topologues, however, is only 2.8. The transitions of CN used 
are blended, with the contributions from the individual hy- 
perfine components impossible to ascertain, thus neither the 
line profiles nor the hyperfine structure ratios could be ex¬ 
amined to constrain the optical depth. It is plausible this 
low ratio is due to the CN transitions being optically thick 
or due to the assumptions made whilst disentangling the 
hyperhne structure not being sufficiently robust. Due to the 
clarity of the line prohles and robustness of the optically thin 
assumption, the HNC and C 2 H ratios are the most reliable 
indicators of ^^C/^®C in IRAS 15194-5115. 


In addition to carbon isotopic ratios in IRAS 15194- 
5115, silicon and sulphur isotopic ratios were determined. 
Due to the high optical depth of ^®SiO, only limits could 
be placed on the ^®Si/^®Si and ^®Si/^°Si. However, both 
^®SiO and ®°SiO appear to be optically thin. The r esult¬ 
ing ^^Si/^°Si of 1.5 is exactly in agreement with Solar. Peng 


et al. (20131 measures the ^®Si/^®Si in 15 oxygen-rich AGB 


stars and finds the sample have values that are either ~ 1.5 
or lower. Few carbon-rich ^®Si/^°Si ratios are reported in 
the literature, due to the low intensity of the ^®SiO and 
^°SiO emission lines. The ^®Si/®°Si has been reported by 


Cernicharo et al. (20001 to be 1.45 ± 0.13 and by He et al. 


I 2008P to helAG ± 0.11 in IRAS 15194-5115. It is cur¬ 
rently uncertain as to whether Si isotopic ratios are affected 
throughout the life of an AGB star. Discussions on this topic 
can be found in Zinner et al. (2006 1 and Decin et al. (20101. 

32s/34s were calculated for both IRAS 15194- 

5115 and IRAS 15082-4808 using two methods: directly from 
CS and C^'^S and indirectly using the ^^C/^^C already calcu¬ 
lated along with ^®CS and C^'^S. In IRAS 15082-4808, both 
methods yield a value of 22, the Solar value. The results 
from IRAS 15194-5115 from the indirect method are almost 
50% greater than that found using the direct method. This 
suggests that the CS line is optically thick - which is also 
suggested by the line profile. As both ®^S and ^"^S are pri¬ 
marily produced by high mass stars in supernovae and the 
32s/34s remained approximately constant at the Solar 
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ratio for the last 7 Gyr (Hughes et al. 20081, this almost- 
Solar ratio in both AGB stars is expected. 


star could be an example of the class of object from which 
the A+B grains originate. 


6.5 Evidence for mixing processes 


Karakas (20101 has published updated stellar yields from 


evolutionary models. The results include the average mass 
fractions of various species in the stellar wind which, as is 
reported, are suitable for comparison with planetary neb¬ 
ula abundances. These models use metallicities of Z=0.02 
(Milky Way metallicity), Z=0.008 (Large Magellanic Gloud 
metallicity), Z=0.004 (Small Magellanic Cloud metallicity) 
and Z=0.0001. The ratio varies between 4 and 

23000 depending upon initial mass and metallicity. At Solar 
metallicity, ratios of less than ten are only found in 

the higher mass progenitors: Matar ^ 5 M©. 

NGC 6537 has been reported in the literature as having 
a core mass of 0.7-0.9 M© and a probable progenitor mass 


of 3-7 Mq (Matsuura et al. 20051. Its C/O ratio has been 
found to be 0.95 ( |Pottasch et al. 2000| . The 5 M© model 
of Karakas ( |2010 1 with Z=0.02 is oxygen rich (C/O = 0.7) 
and has a ratio of 6. The final mass of the star 

is 0.879 Mq. These parameters have reasonable agreement 
with the measured values of NGC 6537, both in this work 
and the literature (e.g. Edwards fc Ziurys|2013 1. This model 
includes hot bottom burning (HBB) but no partial mixing 
zone (PMZ). This suggests that in NGC 6537, hot bottom 
burning was present, in agreement with the findings of |Ed-| 


wards & Ziurys (20131. 


The carbon isotopic ratios of IC 4406 are less reliable 
than that of NGC 6537 due to the low signal-to-noise r atio 
of the ^®C species transitions; Josselin & Bachiller {20031 re¬ 


port '^C=23. The post-AGB mass has been estimated 


by Sahai et al. (19911 to be between 0.6 and 0.76 M©, re¬ 


quiring a lower mass progenitor than for NGC 6537. The 
low-mass progenitor models (1.5-2 M©) of |Karakas| ( |2010[ ) 
at Z=0.02 result in ^^C/^®C of 22 with no HBB or PMZ 
present and a final mass of 0.6-0.64 M©. This implies our 
carbon isotopic ratio is likely a lower limit and that the star 
did not undergo HBB. 


6.6 J-type stars and A-|-B grains 

J type stars, those characterised by ^^C/^^C < 10, have 
been suggested as a potential origin of the SiC presolar 
A-l-B grains ( Zinner et al.|2006 |. The A-l-B grains have low 
i2c/i3c^ a trend between the Si isotopes (known as the 
Si MS line on a plot of d^^ Si against (5^°Si, with equation: 
5^®Si= —20 -I- 1.37 X (5®°Si, Hoppe et al. 201C ’ I and more 


varied S isotopic ratios, but still centred on Solar. 

IRAS 15194-5115 has been measured to have a low car¬ 
bon isotopic ratio, of ~ 4 from column densities of HNC and 
C 2 H, in line with the values obtained from A-l-B grains. The 
silicon ratios are difficult to assess due to the ^®SiO transi¬ 
tion being optically thick. The ®^S/®'*S ratio is also in line 
with that found in A-l-B grains ( [Hoppe et al.|[2Mo| jZinnerj 
et al.|2006 1. This evidence suggests that this probable J-type 


5®Si = 1000 (—') I Decin et al, 
V Si0/^=si0 j 


2010 I 


6.7 Uncertainty considerations 

The uncertainties listed in Tables and 0 are the formal un¬ 
certainties found from propagating the calibration and mea¬ 
surement uncertainties through the calculations. The un¬ 
certainties on the excitation temperature have not been in¬ 
cluded in these calculations as these had to be assumed for 
the majority of the molecular species. In some cases these 
were taken as the closest temperature for which a published 
partition function exists to the temperature resulting from 
the population diagram analyses. Looking at the rotation 
temperatures as a whole, two outliers exists (543 K and 240 
K) and the remainder lie within the range 5-77 K. Taking 
GO as an example: the calculated column density deviates 
from the value at 37.5 K by factors of 2.5 and 3.5 when 
the temperature is changed to 5 K and 150 K respectively. 
Across all molecules, the Scr uncertainty on the calculated 
column densities due to the temperature can generally be 
taken as a factor of 3. 

A number of assumptions have been made in the anal¬ 
ysis within this paper and these assumptions dominate the 
uncertainties of the results. Some of these may be quantified, 
others may not. The hrst is the assumption of optical depth. 
This has been tested by examining the shape of the line pro¬ 
files or, when sufficiently resolved hyperfine structure exists, 
the lines have been fitted to accurately constrain this. Any 
value obtained from a transition deemed likely to be opti¬ 
cally thick is listed as a lower limit to the column density. 

The assumption that the source fills the beam can also 
be tested by examining the line profiles. Some lines show 
resolved profiles and others show unresolved profiles, indi¬ 
cating that the emission regions in each case should not 
deviate too far from a beam filling factor of 1. Therefore 
the assumption that the beam filling factor is unity should 
not contribute significantly to the uncertainties on the cal¬ 
culated column densities. 

The remaining assumptions are difficult to quantify, but 
we estimate that overall the column densities are accurate 
to a factor of 6. 

The isotopic ratios, though dependent upon the column 
densities, are expected to have a lower uncertainty as some of 
the effects discussed above are systematic across the isotopo- 
logues. For example, changes in the excitation temperature 
affect the isotopologues in a similar manner with little effect 
on the ratio of the two. A significant source of uncertainty 
emanates from optical depth assumptions: if the major iso- 
topologue is optically thick then the measured isotopic ratio 
is a lower limit. 


7 CONCLUSIONS 

Presented here are the results of a molecular line survey, cov¬ 
ering the frequency range 80.5-115.5 GHz, carried out over 
a sample of five targets: three carbon-rich AGB stars and 
two oxygen-rich planetary nebulae. The data was obtained 
from the Mopra Telescope, Australia. 

68 individual transitions were identified in the data, 
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emanating from 27 different molecular species. The hyper- 
fine structures of C 2 H and C'^^CH were fitted according to 
the model of Fuller & Myers 119931 to measure the optical 
depth and the contribution to the integrated intensity of the 
blended lines from the individual hyperfine components. Ro¬ 
tation temperatures were calculated for those species with 
multiple detected transitions using the population diagram 
method (Goldsmith & Langer 19991. Source-averaged col¬ 
umn densities were calculated using both the population di¬ 
agram method and a direct calculation using a single transi¬ 
tion. We have also detailed the method used to calculate the 
column densities of species with resolvable and unresolvable 
hyperfine structure. 

From the calculated column densities, carbon, silicon 
and sulphur isotopic ratios have been determined, where de¬ 
tections allow, for the sample. This includes and 

29gi/30si j-atios of the probable J-type star, IRAS 15194- 
5115. This probable J-type star has been shown to exhibit 
isotopic ratios in-line with those expected for the source of 
the A-l-B presolar grains. 
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Figure Al. All transitions in IRAS 15194-5115. Ordinate axis intensities are in units of corrected antenna temperature, abscissa values 
are in units of km/s, corrected for LSR velocity of the source, taken as -15.2 km/s. CO and ^^CO transitions show the stellar emission 
(solid line) and the interstellar contamination (dotted line). 
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Figure A1 - continued 
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Figure A1 - continued 
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Figure A2. All transitions in IRAS 15082-4808. Ordinate axis intensities are in units of corrected antenna temperature, abscissa values 
are in units of km/s, corrected for LSR velocity of the source, taken as -3.0 km/s | |Groenewegen et al.|20d^ . 
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Figure A2 - continued 
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Figure A3. All transitions in IRAS 07454-7112. Ordinate axis intensities are in units of corrected antenna temperature, abscissa values 
are in units of km/s, corrected for LSR velocity of the source, taken as -39.8 km/s ( [Risacher fc van der Tak|2009[ |. 
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Figure A3 - continued 










Figure A4. All transitions in IC 4406. Ordinate axis intensities are in units of corrected antenna temperature, abscissa values are in 
units of km/s, corrected for LSR velocity of the source, taken as -44.5 km/s | |Sahai et al.||199T] |. 
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Figure A5. All transitions in NGC 6537. Ordinate axis intensities are in units of corrected antenna temperature, abscissa values are in 
units of km/s, corrected for LSR velocity of the source, taken as +10.0 km/s < |Edwards Sz Ziurys|201^ . 
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